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LINES OF PRINCIPAL CURVATURE ON CANAL 
SURFACES IN M 3 

RONALDO GARCIA, JAUME LLIBRE AND JORGE SOTOMAYOR 



Abstract. In this paper are determined the principal curvatures and 
principal curvature lines on canal surfaces which are the envelopes of 
families of spheres with variable radius and centers moving along a closed 

rn ■ regular curve in R 3 . By means of a connection of the differential equa- 

tions for these curvature lines and real Riccati equations, it is established 
that canal surfaces have at most two isolated periodic principal lines. 

rl^ Examples of canal surfaces with two simple and one double periodic 

principal lines are given. 

1. Introduction 

The study of principal curvature lines, along which a surface in M 3 bends 
T-i- ' extremely and their umbilic singularities was founded by Monge, Dupin and 

\& • Darboux. See (Gray 1998) and (Sotomayor 2003) for references. For the 

basic facts about principal curvature lines on surfaces the reader is addressed 
to (do Carmo 1976), (Spivak 1979) and (Struik 1988). 

As a consequence of the work of Monge and Dupin the lines of curvature 

^> ■ on quadrics and toroidal -Dupin Cyclides- surfaces were determined. See 

^3 . (Fischer 1986, Chap. 3), for an outline of the theory and for a collection of 

remarkable illustrations. 

In (Gutierrez and Sotomayor 1982, 1991, 1998) ideas originating in the 
Qualitative Theory of Differential Equations and Dynamical Systems, such 
S^ ■ as Structural Stability and Genericity, were incorporated into the subject; 

see also (Garcia and Sotomayor 2002). Historical comments, going up to 
some recent developments on principal curvature lines, can be found in (So- 
tomayor 2003). 

The global dynamic complexity of principal curvature in simple smooth 
surfaces was illustrated in (Gutierrez and Sotomayor 1991, 1998). This in- 
cludes examples of recurrent (dense) principal curvature lines on spheroidal 
surfaces, which are small perturbations of ellipsoids of revolution and canal 
surfaces with constant radial functions. The methods established in these 
works show how to make hyperbolic a principal cycle by means of a small 
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smooth perturbation. This leads to smooth toroidal immersions with arbi- 
trary large number of principal isolated principal cycles. 

In this paper is improved and reproved a result that goes back to Vessiot, 
establishing that principal curvature lines on canal surfaces immersed in M 3 
(see Section |21 Definition EJ) verify a Riccati equation (see Sectional Remark 
EJ). In fact, (Vessiot 1919) is the first source known by the authors for the 
connection between Riccati equations and principal curvature lines on canal 
surfaces. The improvement consists in the formulation of precise conditions 
for canal surfaces be regular immersions and also to have umbilic points 
(see Section El Theorem |3] and Remark EJ) . A consequence of the Riccati 
structure for principal curvature lines on canal immersed surfaces implies 
that the maximal number of isolated periodic principal lines is 2. Examples 
of canal surfaces with two (simple i.e. hyperbolic) and one (double i.e. 
semi-stable) principal periodic lines, absent in (Vessiot 1919) and also in 
later references, are given in this paper (see Section^ Proposition IHJ. 

2. Regular Canal Surfaces 

Consider the space M 3 endowed with the Euclidean inner product < , > 
and norm | | =< , > 1 / 2 as well as with a canonical orientation. The wedge 
product A of vectors is defined relative to this orientation. 

Let c be a smooth regular closed curve immersed in M 3 , parametrized by 
arc length s € [0, L]. This means that 

c'(s) = t(s), \t(s)\ = l, c(L) = c(0). (1) 

Assume also that the curve is bi-regular. That is: 

k(s) = \t'(s)\ > 0. (2) 

Along c is defined its moving Frenet frame {t,n, b}. Following (Spivak 
1979) and (Struik 1988), this frame is positive, orthonormal and verifies 
Frenet equations: 

t'(«) = n(s)n(s), n'(s) = -«(a)t(«) + r(s)b(s), b'(s) = -r(s)n(s). (3) 

Equations (fT ]) to P | ) define the unit tangent, t, principal normal, n, cur- 
vature, k, binomial, b = t A n, and torsion, r, of the immersed curve c. 

Proposition 1. Letr(s) > and8(s) € ]0, n[ be smooth functions of period 
L. The mapping a : T 2 = S 1 x S 1 -* K 3 , defined on R 2 modulo L x 2n by 

a(s, (p) = c(s) +r(s) cos#(s)t(,s) + r(s) s'm9(s)[cos(pn(s) + sin</?b(s)], (4) 

is tangent to the sphere of center c(s) and radius r(s) if and only if 

cos0(s) = -r'(s). (5) 

Assuming ©, with r'(s) < 1, a is an immersion provided 

«(,) < '-r'M' -'(.»*(.) (6) 

r(s)s/l — r'(s) 2 
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where 



Proof. Calculations using equations © give 

a s = — - = alt(s) + a 2 n(s) + agb(s), 
as 

a s = 1 + r' cos 9 — r9' sin 9 — rn sin 9 cos ip, 

a s = r' sin 9 cos tp + rO' cos cos tp — rr sin 9 sin 93 + rs cos 9, (7) 

a s = r' sin sin 93 + r9' cos sin ip + rr sin cos 99. 
Here, of course, r, 0, r, k are functions of s. Also, 

a v = — — = r{s) sin#(s)[— sin</?n(s) + cos^b(s)]. 

The unit normal vector pointing inward the sphere \p — c(s)\ = r(s) with 
center c(s) and radius r(s) at p = a(s,<p) is N a = (c(s) — a(s,<p))/r(s), 
which is given by: 

N a = — cos 9(s) t(s) — sin0(s)[cos</3n(s) + sin<£>b(s)]. (8) 

Clearly, < a^^Na >= 0. Calculation gives < a s ,N a >= — (cos# + r'). 
Therefore the condition of tangency at (s, tp) of a to the sphere \p — c(s)| = 
r(s) is cos6>(,s) = — r'(s). 

Additional calculation gives: 

F(s, ip) = < a s ,a v >= rr 2 sin 2 9 + K(s)r'r 2 sin 9 sin 99, 

(9) 
G(s, ip) =< a^, a v >= r 2 sin 2 = r 2 (l — r' ). 

At this point it is appropriate to write the mapping a and its derivatives 
involving only functions of r(s) rather than of 9(s), replacing the expressions 
cos 9(s) = -r'(s), sin0(s) = (1 - (r') 2 ) 1 / 2 in equations © and 0. 

The expression for E is as follows: 

E(s, ip) = < a s ,a s >= k r (1 — r' ) cos (p + 



2 nr[rr' — (1 — r' )/v 1 — r' ] cos 99 + 

2nTr 2 r'Vl - r' 2 sin tp + (1 - r' 2 ) + 2rr"+ ^ 10 ^ 

f r ">'- ' 2 \2 1 „/2. 2/i _/2s , „//2! 



-g-T (1-r ) +r k (1-r )+r , 

(1 — r ) 

Additional calculation and simplification using equations Q, with the 
tangency condition Q imposed, and (jlT)|) gives 



EG-F 2 = r{s) A yi-r'{s) 2 K{s) cos <p + r"(s) - (1 - r'(s) 2 )/r(s)} 2 . (11) 

This expression, which is equal to \da/ds Ada/dtp\ 2 , vanishes if and only if 

1 - r'{s) 2 - r{s)r"(s) 

cos<p = — = . (12J 

y 7 ! — r'(s) 2 K(s)r(s) 
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Condition © states that the absolute value of the right-hand member of 
l)12j) is larger than 1, which by IjllJI implies the linear independence of a s 
and a» at every point (s, (p), namely that a is an immersion. □ 

Definition 2. A mapping such as a, of T 2 into M 3 ; satisfying conditions 
@ ond (jHJ) m/Z be called an immersed canal surface with center along c(s) 
and radial function r(s). When r is constant, it is called an immersed tube. 
Due to the tangency condition l(5|). the immersed canal surface a is the 
envelope of the family of spheres of radius r(s) whose centers range along 
the curve c(s). 

Remark 3. Proposition^ is partially found in (Vessiot 1919). There, how- 
ever, the regularity condition © was overlooked. Other references for canal 
surfaces are (Blaschke 1929) and (Gray 1998). 

3. Principal Lines on Regular Canal Surfaces 

In this work the positive orientation on the torus T 2 = S 1 x S 1 is defined 
by the ordered tangent frame {d/ds,d/dip}. Therefore the positive unit 
normal -or Gaussian map- N a of the immersion a is defined by \da/ds A 
da/dip\N a = da/dsAda/dip. By the tangency condition ® this unit vector 
is given by (|%|). which can be written as follows: 

N a = r'(s)t(s) - (1 - (r'(s)) 2 ) 1/2 [cos93n(s) + sin ¥ ?b(s)]. (13) 

It points inwards the toroidal surface defined by a. 

Below will be studied the global behavior of the principal curvature lines 
of a. 

Theorem 4. Let a : S 1 x S 1 — > IR 3 be a smooth immersion expressed by Q . 
Assume the regularity conditions © and © as in Proposition^ and also 
that 

l-r\sf-2r(s)r"{s) 

2r(s)y/l-r'(s) 2 ' 

The maximal principal curvature lines are the circles tangent to d/dip. The 
maximal principal curvature is 

k 2 {s) = l/r{s). 

The minimal principal curvature lines are the curves tangent to 

d_ 

ds~ y y "' ' (] 

The expression 

r(s)(l - r'(s) 2 ) 1/2 K(s) cos tp + r(s)r"(s) - (1 - r'(s) 2 ) (16) 



V(s, ¥>) = £-( r(s) + _ ^ } )2)1/2 k(s) sin y) ±. (15) 
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is negative, and the minimal principal curvature is given by 

kl ( s(p ) = k(s)(1 -r>(s) 2 )V 2 cos <p + r"(s) 

r(s)(l — r'(s) 2 ) 1 / 2 K(s) cos ip + r(s)r"(s) — (1 — r'(s) 2 ) 
There are no umbilic points for a: fci(s, (/?) < /^(s). 

Proof. Direct calculation gives: 

<9iV Q 5a 

By Rodrigues equation, see (do Carmo 1976) and (Struik 1988), the circles 
s = constant are principal curvature lines of a, with principal curvature 
k 2 {s). 

Denoting the metric of a by Eds 2 + 2Fdsdip + Gd(p 2 , it follows that the 
direction orthogonal to d/dip, giving the other principal direction, is defined 
by the vector field 

ds dtp' 
which by equation Q is collinear with V(s,tp) in (|15j) . 

Differentiation of (fT3|) using Frenet equations (JHJ) gives 



and 



W* = y + K (i - r ' 2 )V2) cos ]t(s)+ 

OS 

[(r'r" cos ip + r(l — r' ) sin ip-\- 

r'«(l - r' 2 ) 1/2 )/(l - r ,2 ) 1/2 ]n(s)+ 
[(-r(l-r' 2 ) 1 / 2 cos99 + rV / sin(^)/(l-r' 2 ) 1 / 2 ]b(s), 



(18) 



\J\ — r'(s) 2 [sin ip n(s) — cos </? b(s)] . 
Substitution, taking into account equation (|15|) . leads to 



dN a . dN a 

-W {S ^ ) = ^s- 



T+ {l _ r/ 2 ) l/2 KS[n ^ 

Similarly for a, calculation of its derivative ^fe, gives 



dNa 

dip 



(19) 



(jcy. t (s) 

—(s,(p)=A(s,(p)[t(s) + - 2 (cos ipn(s) + sinyb(s))], 

c»l/ (1 - r'(s)) 1 / 2 (20) 

^4(s, u) =1 — r'(s) — r(s)r"(s) — K(s)r(s) cos p(l — r'(s) ) ' 

Further calculation using equations (|18jl , (|19() , Q2U|) and the expression for 
A;i(s,(/?) given in equation ((T7)) . follows that 

dN a . , , , ,9a, . 
-^ r (s,(^) + A;i(s,v7) — (s,(^) = 0. 
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This is Rodrigues equation, which establishes that, in fact, k\ is a principal 
curvature. 

From © it follows that (|16ft . the denominator of fci, is always negative. 
This implies that k\ < l/r(s) = /c2- Otherwise, k\ > &2 and, after direct 
manipulation, by ([Fiji this would lead to 



cos If 



> 



1 - r'(s) 2 - 2r(s)r"(s) 



2«;(s)r(s) v / l - r'(t 



> 1. 



(21) 



This prevents the existence of umbilic points and justifies the names maxi- 
mal, for subscript 2, and minimal, for subscript 1, given in the statement. □ 

Remark 5 (Umbilic Points in Canal Immersions). The calculation ofk\ and 
the condition for the appearance of umbilic points for regular canal surfaces 
has not been considered in previous works on the subject. The discussion 
leading to ()21|) leads to the following equation for umbilic points: 



1 



COS(/2 



r'(s) 2 



2r(s)r"(s) 



2K(s)r(s) v / l-r / ( 



It gives a non-empty curve if condition (|14|) is not imposed. This curve 
however consists of removable singularities for the principal line fields which 
have smooth extensions to the whole torus, given by d/dip and V. Under 
generic conditions on the radial function, r(s), and curvature, k(s), func- 
tions it will be expected to appear curves of umbilic points as those studied 
in Proposition 2 of (Garcia and Sotomayor 2005). 

Remark 6 (Minimal Principal Foliation in terms of Differential Forms). 
Written as a differential form, the vector field ()15|) becomes: 



uj\ = dip + 



t(s) + 



r'(s) 



(1-^(5)2)1/2 



n(s) simp 



ds. 



(22) 



When it is used another parameter t, say T-periodic, related with the arc 
length s of c in the form s = s(t), with s =< c,c > ' , the form in (|22j) 
writes as: 



oj\ = dip + 



r(t) + 



r{t) 



(s(t) 2 - r(t) 2 )V2 



n{t) sinyj 



s(t)dt. 



(23) 



Remark 7 (Vessiot). By the change of coordinates tan((^/2) = z, equation 
((2*2*|) is transformed into the L-periodic Riccati equation 



1 



t(s)(1 + z 2 ) + cot 6(s)k(s)z. 



(24) 



Therefore, the solutions of equation (|15|) can be obtained from the solutions 
of equation ()24|) contained in [0,T] x [— ir, it]. 
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4. Canal Surfaces with one and two Principal Cycles 

In this section is carried out a discussion on the qualitative properties of 
equation (|15|). or of its equivalent form (|23j) . absent in (Vessiot 1919). 

Being equivalent to a periodic Riccati equation, this equation has a Moebius 
transformation as return map (see Hille 1976) and therefore can have either: 

(a) all its solutions periodic, 

(b) all its solutions dense, 

(c) two hyperbolic (simple) periodic solutions or 

(d) one semi-hyperbolic (double or semi-stable) periodic solution. 

An example of situation (a) is exhibited by the standard torus of revolu- 
tion. An example of (b) is given in (Gutierrez and Sotomayor 1991, 1998), 
for a canal surface of constant radius (a tube) around a curve that is not 
bi-regular. Below will be given examples of cases (c) and (d). This will also 
provide examples of case (b) for bi-regular curves. The cost of this is heav- 
ier calculation which, nevertheless, is easy to corroborate with Computer 
Algebra. 

The example consists in a deformation of the T = 2-n-periodic plane 
elliptic curve 

c(t) = (2 cost, sint, 0), 

whose curvature k(i) = n(t,0) is k(£) = 2/(4 — 3cos 2 t) 3 / 2 . 

Proposition 8. Consider the three parameter family of canal surfaces S eyPyfJi 
around the curve c £ (t) = (2 cost, sint, ek{t)), with radial function 

r(t,/i) = p + pk(t), k(t) = -18 cost sint/ (4 - 3 cos 2 t) 5/2 . 

There are two smooth curves e = £\(p) = —(/, + 0\{p 2 ) and e = £2^) = 
/i + 02(/U 2 ) ; such that for any p, p small and positive, the canal surface S £jPjfl 
has two hyperbolic principal cycles for e e]ei(/z), £2 (/•*)[ and has one double 
principal cycle along the curves e = £\{p) and e = £2(1^)- 

Proof. Standard calculation gives: 

s{t,e) = (4sin 2 t + cos 2 t + (£K(t)) 2 ) 1/2 = (-3 cos 2 t + 4) 1 / 2 + 0(e 2 ), 

for the element of arc length of c e . 
Also, the curvature of c £ is 

' 1458 cos 12 t + 486 cos 10 t - 4671 cos 8 1 



K(t,e) = K(t,0) + 648£ 



(4-3cos 2 t) 19 / 2 
432 cos 6 t + 4296 cos 4 1 - 2016 cos 2 t + 16" 



+ 



(4-3cos 2 t) 1 9/ 2 +0{e). 

Here will be needed only k(t, 0) = k{t) as given above and k{t) in the integral 

2n , x ., . , /' 27r 36(9cos 4 t-4cos 2 t-4) 8829 . . 

s(t,0)K(t)n(t)dt = j o _L_____^ (ft = __„. (25 ) 
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The torsion of c £ is 

dr 
r(t, e) = r(t, 0) + s— (t, 0) + 0(e 2 ). 

Here will only be needed the expressions for r(t, 0) = and T\(t) = dr/de(t, 0). 
From direct calculations follows that 

108 (54 cos 8 t + 207 cos 6 t - 369 cos 4 t + 68 cos 2 t + 44) 
Tl(t) " (4-3cos 2 t) 11 /2 ' 

f 2n 8829 

C = -J i(t J 0)7i(t)A = ^ S 7r. (26) 

The differential equation corresponding to (|23|) for the 27r-periodic canal 
surface S £<p ^ is 

dg9 TTr . , . / ak{t)n(t, e)sm.(p \ , _. 

I = *ft *>. -. rt - .<*, .) (-(*. .) - ^rkw) ■ (27) 

Denote by Tl(ipo,£, fi) the return map for equation (|27|h It is given by 
II(</?0) £) A*) = $(271-, (^o, e, At), where $(£, po, e, A*) is the solution of (|27|) such 
that <J>(0, ipo,e,/j,) = tpo. Therefore, the 27r-periodic solutions are given 
by the implicit surface Ii(ifo,£, /u) — ipo = 0. Since n(<^o,0,0) = ipo, the 
Fundamental Theorem of Calculus (Hadamard's Formula) implies that 

f 1 dU f 1 dn 

U((p ,e,n) -cp = £ / —(tpo,u£,uiJ,)du + ii -—(ip G ,u£,un)du. (28) 

Jo oe J dn 

Write 

on <9n 

— (ipo,u£,ufi) = —{ipo,0,0) + P{ip ,u,£,fi)£ + Q(ip ,u,£,n)i2,(29) 

dU on 

— (ipo,ue,ufi) = —(ip ,0,0) + R(ip ,u,e,fi)e + S(cp ,u,£,n)iJ 1 .(30) 

The expressions for the derivatives of the solutions of the differential 
equations with respect to parameters applied to (|27|). which in the present 
case, following classical differential equations results, are the integrals of 
the non-homogeneous linear — called variational — equations. According 
to (Sotomayor 1979, page 42), (Chicone 1999, page 337) or (Coddington 
and Levinson 1955, page 30), where these classical results are proved, these 
equations and their initial conditions are as follows: 

® e (t,<po,0,0) = W v (t,<po,0,0)* B (t,<PoA0)+W e (t, <po,0,0) 

= -s(t,0)n(t), $ e (0,v>o,0,0) = 0, 

%(t,(po,0,0) = ^(^o,0,0)^(t,^o,0,0) +Wyt, pc 0,0) 

= -s{t,0)k{t)K(t)sinip , <^(0, (p Q ,0,0) = 0. 
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Integrating these equations, taking into account equations (|23)l and Q26|>. 
leads to 

■(^o,0,0) = $ £ (2vr,^ ,0,0) = - / ri(v)cfo = C, (31) 

JO 



<9e 



dll 

— (99o,0,0) = ^(2vr,^ ,0,0)=C7sin^ . (32) 

Writing e = v\i and substituting into (|28|). taking into consideration the 
expressions (|3~Tj) and (|32() . obtain 

n((^o, e, m) - <£o = C/i(z; + sin (p + i/Z(i/, /U, (p )). 

Therefore, the surface of 27r-periodic orbits, in the (y , /j,, (£>)-space consists 
on a plane \i = crossing transversally a regular sheet P at the sinusoidal 
curve of equation v + sin ipo = 0. The critical values of the projection on the 
(i>, fi)— plane of the surface P, coming from double (semi-stable) periodic 
solutions, the fold curve, must cross transversally the z^-axis at points cor- 
responding to v = ±1, which are the critical values of the projection of the 
sinusoidal curve. Going back to the coordinates (e, /u) gives the conclusion 
formulated in the proposition. □ 

Remark 9 (Arnold Tongues). Due to equations (|31jl and (|27jl. for p small, 
the return map on the tube S EiPt o is a rotation with its rotation number 
changing monotonically with e (taken small). Therefore it takes irrational 
values and the tube presents case (b). It also takes rational values, q/p, which 
are the vertices of the Arnold Tongues (in the (e , jj) -plane) , corresponding 
to canal surfaces with periodic closed principal lines winding p-times around 
the parallels (t-circles) and q-times around the meridians (cp-circles) of the 
torus. The sector with vertex at (0, 0) established in Proposition [3 is the 
principal tongue, with rotation number 0. See (Chicone 1999, page 372). 
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